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Abstract 

Recently it has been found that the structure of Skyrmions has a close analogy to 
that of fullerene shells in carbon chemistry. In this letter we show that this analogy 
continues further, by presenting a Skyrme field that describes a lattice of Skyrmions 
with hexagonal symmetry. This configuration, a novel 'domain wall' in the Skyrme 
model, has low energy per baryon (about 6% above the Faddeev-Bogomolny bound) 
and in many ways is analogous to graphite. By comparison to the energy per baryon 
of other known Skyrmions and also the Skyrme crystal, we discuss the possibility of 
finding Skyrmion shells of higher charge. 
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1 Introduction 



Our recent work [0], involving numerical simulations of the full nonlinear field equations 
of the Skyrme model, has revealed that the structure of minimal energy multi-Skyrmions 
has a rich and fascinating complexity, where surfaces of constant baryon density are given 
by trivalent polyhedra with holes at the centre of each face. These structures are very 
reminiscent of those appearing in fullerene chemistry to describe closed shells of carbon 
atoms P, ^, leading us to draw an analogy between Skyrmions and carbon chemistry. 

More explicitly, the empirical Geometric Energy Minimization rule of ref. [0] states 
that, the baryon density isosurface of a charge B Skyrmion is an almost spherical trivalent 
polyhedron with 4{B — 2) vertices, 2{B — 1) faces and Q{B — 2) edges. For B > 7 the 
polyhedron comprises twelve pentagons and 2{B — 7) hexagons, which is precisely the 
structure of a fullerene corresponding to a closed shell containing 4(5 — 2) carbon atoms. 
In particular, the B = 7,8 and 9 minimum energy configurations were seen to have the 
same structure as shell-like forms of C20, C24 and C28 respectively fl], |]. Motivated by this 
numerical work, a new approach to constructing Skyrme fields 0, based upon rational 
maps between Riemann spheres was developed. This has allowed a good mathematical 
understanding of some aspects of these Skyrmions and has also produced a Skyrme field 
of charge seventeen which has the structure of the most famous of the fuUerenes, the Cgo 
Buckminsterfullerene. It too appears to be the minimal energy Skyrmion of this charge. 

In considering very large fuUerenes, where hexagons are dominant, the twelve pentagons 
may be viewed as defects, inserted into a flat hexagonal structure, in order to generate the 
required curvature necessary to close the shell. Energetically the optimum structure is an 
infinite hexagonal lattice, that is, a sheet of graphite; the most stable form of elemental 
carbon from the thermodynamic point of view. The reason that closed shells are preferred 
for a finite number of carbon atoms is that the penalty for introducing the pentagonal 
defects is not as severe as that incurred by having dangling bonds at the edges of a truncated 
graphite sheet. A prediction of the fullerene approach to Skyrmions is therefore that a 
Skyrme field should exist which represents a hexagonal lattice, that is, the analogue of 
graphite. Furthermore, although this configuration would have infinite energy, since it has 
infinite extent in two directions, its energy per baryon should be lower than that of any of 
the known finite energy Skyrmions. In the next section we shall verify this prediction and 
thus provide another piece of evidence in support of the similarities between fuUerenes and 
Skyrmions. 

For completeness, we include expressions for the Lagrangian of the Skyrme model, 
which in terms of the su(2)-valued right currents i?^ = {d^U)U'^ is 

C = -iTr(i?^i?'^) - ^Tt{[R^,R,][R>^,R'']), (1.1) 

where we have used scaled units of energy and length, and the baryon density B, whose 
spatial integral gives the integer-valued baryon number B, is given by 

B = --^e,,kTiiR,R,Rk) . (1.2) 
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As throughout this letter, we take latin indices to run over the spatial values 1, 2, 3. The 
corresponding Faddeev-Bogomolny bound on the energy E is simply E > \B\ x 127r^. 

2 An ansatz with hexagonal symmetry 

Following a suggestion of Atiyah, it has recently been proved by Jarvis 0] that the space of 
SU(2) BPS monopoles is diffeomorphic to the space of equivalence classes of rational maps 
between Riemann spheres. These rational maps arise from the monopole as the scattering 
data of a linear operator when considered along all possible lines emanating from a chosen 
origin. Through the use of a new ansatz for Skyrme fields 0, it has been possible to 
use Jarvis rational maps to construct good approximations to the known minimum energy 
Skyrmions. The fact that Jarvis maps, being the scattering data along radial lines, are 
relevant reflects the feature that these Skyrmions are shell-like structures. 

In this letter we are concerned with a two-dimensional lattice of Skyrmions, which is 
clearly not a shell-like configuration, making Jarvis maps inappropriate for this application. 
However, an older diffeomorphism of Donaldson [^], between SU(2) monopoles and based 
rational maps, is appropriate. These Donaldson maps arise as the scattering data of a 
linear operator considered along all lines in a chosen direction. This definition requires a 
decomposition of into R x (D, which is exactly the situation for a lattice configuration, 
where the complex plane refers to the plane of the lattice and the real coordinate is the 
height above the lattice. To be precise, since we wish to consider an infinite lattice, the 
Donaldson map will not be rational, but instead is required to be merely meromorphic; 
which may be regarded as an infinite limit of a rational map. 

A consideration of the ansatz introduced in ref. 0, together with the modifications 
discussed above, leads us to the following Skyrme field ansatz 

U{x,, X2,x,) = exp (^.^-^L_{Wt. + Wt+ + (1 - \W\^)n)^ , (2.1) 

where denote the Pauli matrices with t± = ti ±ir2, W G(DP^ is a holomorphic function 
of z = xi + ix2, and / G R is a function of x^. 

The lattice occupies the xiX2 plane, in which we use the complex coordinate z. Thus 
from the above ansatz we see that the direction of the vector of pion fields is determined 
by the (DP^ field W, given the position in the lattice, whereas the length of the vector of 
pion fields is determined by the profile function /, given the height above the lattice. 

The next issue to address is that of the boundary conditions on the Skyrme field U 
and hence on the functions / and W. To have a periodic lattice there must exist complex 
constants Qi and ^2 which are the fundamental periods of the lattice, that is, 

U{z + 71^1 +mn2,X3) = U{z,X3) Mn.m^TL. (2.2) 

Let denote the associated fundamental parallelogram, that is, the torus given by the 
region in the complex plane with vertices 0, i^i, ^2, + ^2 and opposite edges identified. 
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We can now restrict our analysis to the region of given by R x T^, with the field in 
the remaining regions determined by periodicity. Thus we see from our ansatz that 
W is required to be a holomorphic map : h-> (DP . 

To consider the boundary conditions in the direction orthogonal to the plane of the 
lattice we need to recall our motivation. The lattice is being thought of as an infinite limit 
of the shell-like Skyrmions containing pentagons and hexagons. Thus, in approaching this 
limit, we imagine the lattice as being a part of the bottom of a larger and larger shell 
and hence below the lattice is the outside of the shell, where the Skyrme field tends to 
the vacuum, giving the boundary condition lim2,3^_oo U = I2. However, above the lattice 
is the inside of the shell, where the Skyrme field is approaching the negative vacuum 
associated with the centre of the Skyrmion, so the appropriate boundary condition is 
limx^-^+ooU = — 12- Note that this implies that our Skyrme lattice is a domain wall, 
separating regions of differing vacuum values. Examination of our ansatz now reveals 
that the boundary conditions for the profile function /(X3) read 

/(-oo)=0, /(oo) = 7r. (2.3) 

To compute the baryon number and energy of the Skyrme field ( [2.1D in a fundamental 
section we shall follow the approach of ref. 0. The strain tensor, defined as 

A, = ~Tt{RiRj) , (2.4) 

is symmetric and positive semi-definite. If it has eigenvalues A^,A2,A| then the Skyrme 
energy density S and baryon density B are given by 

S =Xl + Xl + Xl + Xlxl + Xlxl + Xlxl (2.5) 
B = AiA2A3/27r2. (2.6) 

For the ansatz (|2.1| ) the strain in the direction normal to the lattice is orthogonal to the two 
strains in the directions of the lattice, which are equal. Therefore the Aj may be interpreted 
as the strains in the Xj directions, making it is easy to show that 

A3 = /', Ai = A2 = 2Jsin/, (2.7) 

where we have defined the quantity 



J = 







1 + 


\W 


2 



(2.8) 



Substituting the expressions for the strains ( ^.T] ) into those for the energy and baryon 
density, we arrive at the result 

S =/'2 + 8J'(/'2 + sinV) + 16J'sinV, (2.9) 
B =AjYsinV. (2.10) 
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We now wish to use these densities to compute the energy E and baryon number B 
in a fundamental section of the lattice, by integrating over the region X3 G (—00, 00) and 
{xi,X2) e T^. To do this we note that since IV is a map W : T"^ (DP^, then it has 
an associated integer, k, which is its degree. Explicitly, k is given by integrating over the 
torus the puUback under W of the Fubini-Study area 2-form on (DP^, which in this case 
gives 

k=-f J2 dxidx2, (2.11) 

since is a holomorphic function of z. 

Using (|2.10|) it is now easy to see that the baryon number is equal to the degree k, since 



2 r 
B = — f sin^ / (ix3 / dxidx2 

TX^ J -00 JT^ 



/-^sin2/ 



k, (2.12) 



where we have used the expression ( p. 11 ) and the boundary conditions (|2.3| ). 



In calculating the energy it will be useful to introduce a scale parameter /i by writing 
u = X3/ and setting /(xa) = g{u). Then, if A is the area of the fundamental torus T^, 
integrating the density (Of) gives 



where the EiS are the following integrals over u 



/oo r A 1 U 

dx3 / dxidx2 £ = -El + -E2 + ^lE^ + ^E^ (2.13) 
-00 JT2 an A 



El 


J — 00 


g'^ du , 

roo 


E2 


= Silk 


/ 

J —CO 

roo 


E3 


= Sirk 


/ sin 

J ~oo 
roo 


E4 


= 16J 


/ sin 

J —00 



(2.14) 



The only remaining dependence on the map W is the quantity I, which is defined as 

T = A I dxidx2, (2.15) 

and has the important property that it is independent of A. 

The scale fi and area A can now be determined, in terms of the -Ej's, by minimization 



of the energy (|2.13| ). Requiring = ff = 0, gives the result 



/i = ^E2/E3 , A = ^E2EjEiE3 , (2.16) 
and hence a minimized energy of 



E = 2{^EiEi+ ^E2E3). (2.17) 
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To proceed further we now need an explicit expression for the map W{z). To obtain a 
holomorphic map from the torus, we take W to be an eUiptic function of z. Exactly which 
elliptic function to take is determined by the fact that we wish to construct a hexagonal 
lattice, so we require a fundamental period parallelogram which has a 60° angle between 
the two fundamental periods. The appropriate elliptic function is the Weierstrass function 
p{z) satisfying 

p" = ^p'-l), (2.18) 

which has real period Vli = r(|)r(|)/(2-\/37r) and imaginary period given by Q2 = ^i(l + 
iV3)/2. From this we see that the period lattice is equilaterally triangular, and thus we 
have the desired 60° angle. Obviously we can scale both the elliptic function and its 
argument and still retain the above desired property, hence we take 

W{z) = cp{z/a), (2.19) 

where c and a are arbitrary real constants. Note that by the inclusion of the factor a we 
must now rescale the fundamental torus. Furthermore, for later computational purposes 
it is convenient to work with a rectangular fundamental torus, which is achieved by taking 
it to be the given by {xi,X2) G [0,afii] x [0,av^fii]. As this torus contains two 
fundamental parallelograms and the p-function has a double pole in each of these then, by 
counting preimages, we see that the degree of the map in this case is = 4. 

For a given k, the energy E is minimized by minimizing the value of the integral X 



defined in (|2.15|) . This is clear since its only appearance is as a coefficient in front of the 
positive term in ( |2.14| ). Recall that X is independent of the area A, and hence a since 
the two are simply related by 

A = VSa^nj . (2.20) 

Computing X for the one-parameter family of maps ( |2.19| ), given by varying c for any fixed 
a, we find that its minimum value is I ^ 193, which is attained when c ~ 0.7. 

In order to continue with an analytical treatment we now make an ansatz for the profile 
function g{u), which we shall see turns out to be a reasonably good choice. We choose the 
sine-Gordon kink profile function 

g{u) = 2arctan e"" , (2.21) 



which has the advantage that all the integrals in ( p. 14 ) can be performed exactly. The 
results are 

El = 2 , E2 = 1287r/3 , E3 = GAtt , E^ = 64X/3 , (2.22) 
from which we find that the scale and area are 



fi=\/2/3, A=^^, (2.23) 



and using ( p.l7|) the energy is 



^ = 16A/2/3(VX + 87r) . (2.24) 
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Recalling that B = k = A we thus compute that the energy per baryon is 



E/B = 1.076 X 127rl (2.25) 

In discussing the energy of various configurations it will be useful to define A to be the 
percentage excess over the Faddeev-Bogomolny bound ie. 



E - 127r2 


B\ 


127r2 


B 





A = ' ' X 100%. (2.26) 



Thus, our ansatz for the lattice has A^- = 7.6%. 



3 Skyrmion Architecture 

In this section we discuss the accurately computed energies of various configurations, using 
a numerical relaxation of the full nonlinear Skyrme model. Using this information we 
speculate on the kinds of structures which may form for Skyrmions of high charge. 

In the previous section we were able to compute an explicit Skyrme field which describes 
a lattice with a relatively small excess energy. However, some approximations were made, 
for example, the choice of the profile function (|2.21| ) and so the energy of the true lattice 



will be lower than the ansatz value of A,~ = 7.6%. To determine the true value, we take 

lat ' 

the ansatz field as a starting configuration in a numerical relaxation computation using the 
code described in detail in ref. 0. The simulations were performed on a grid containing 
100 X 100 X 58 points and periodic boundary conditions imposed in appropriate directions. 
The computation of the initial excess confirmed the value to be A = 7.6% but after 
relaxation, which also involved some minor rescaling to be sure to obtain the minimum 
energy lattice, the excess of the true lattice was found to settle down to Aiat = 6.1%. 

In Fig.l we display a surface of constant baryon density for the hexagonal lattice. The 
hexagonal structure is clearly visible, with the baryon density isosurface having a hole in 
the centre of each of the hexagonal faces. Note that the section of lattice we are considering 
contains exactly eight full hexagons and is of baryon number four, computed to be 5 = 3.84 
on the discretized grid. The fact that each hexagon may be thought of as having baryon 
number one-half is the expected infinite limit of the polyhedron structure discussed in the 
introduction, where a charge B Skyrmion is found to have 2(5 — 1) faces. 

As mentioned in the introduction, the excess energy of the lattice is lower than that 
found for any finite charge Skyrmion. For a single Skyrmion the excess is Ai = 23%, 
whereas for the configurations up to charge nine the minimum excess occurs for i? = 9 
which has Ag = 9.8% fl]]. A relaxation of the charge seventeen buckyball configuration 
gives the smallest known value A17 = 7.2% Thus our results are all consistent with 
the fuUerene picture of Skyrmions where, at least for Skyrmions of modest charge, the 
structure is a shell of pentagons and hexagons. 

The fact that Aiat is so low is encouraging for the possibility of Skyrmion shells at 
even higher charge, however it is not low enough to conclude that shell-like structures 



7 



continue indefinitely. This is because of tlie existence of the Skyrme crystal 0, which 
is a configuration that is periodic in all three space dimensions and consists of a crystal of 
half-Skyrmions. A Fourier series analysis approximates the energy excess of the crystal 
to be A~ = 3.8% and a simple analytical formula exists which gives a good approximation 
to the fields of the Skyrme crystal 0. Using this as an initial configuration in our relaxation 
scheme, we find a true energy for the Skyrme crystal of Aery = 3.6%. A surface of constant 
baryon density is shown in Fig. 2. The half-Skyrmion structure of the crystal is evident and 
it is clear that this configuration is of a very different type to that of the lattice shown in 
Fig.l. 

Naively, the fact that Aery < Aiat suggests that there could be some value of the 
Skyrmion charge, -Berit, at which the shell-like structures will be replaced by more three- 
dimensional configurations, but the details of how and when this might take place are 
unclear. However, a simple comparison of Aery and Aiat measures the volume effect, but 
since these two values are not so different the crucial factor will be an area effect, which is 
associated with the fact that a finite portion of the Skyrme crystal needs to be smoothed 
off at the edges. One suitable candidate for completing the edge of the Skyrme crystal is a 
configuration like the face of the B = 4 cubic Skyrmion. However, since this has an excess 
of A4 = 12%, quite a substantial piece of the crystal needs to be included before such a high 
penalty for its boundary could be accommodated. In contrast, a large shell structure only 
ever requires twelve pentagon defects and they can be included at a relatively small cost, 
as demonstrated by the B = 17 buckyball with A17 = 7.2%. All this suggests that even if 
the Skyrme crystal structure appears at some charge -Berit then this may well be very high. 
Therefore, there is at least the possibility of a range of modest charges where other fullerene 
structures may exist. Other exotic possibilities include analogues of the bucky-tubes (long 
thin configurations comprised of spirals of hexagons with caps containing pentagons) and 
also shells inside shells. 

If the above structure change to the Skyrme crystal does indeed occur then for B > Berit 
then the fullerene analogy will be lost, since the Skyrme crystal is a configuration with 
valency six. However, it is possible that the known Skyrme crystal is not the minimum 
energy crystal structure. If, for example, a configuration exists with the structure of 
the diamond lattice, then the similarities with carbon atoms could be maintained. This 
is a four-valent lattice with tetrahedral symmetry, but would not be as simple as the 
tetrahedral lattice formed from individual Skyrmions in the nearest neighbour attractive 
channel, which relaxes to the cubically symmetric Skyrme crystal. However, it is a difficult 
task to investigate this possibility since a Skyrme field with the correct properties needs 
to be found before numerical simulations can be performed. Using the ideas from rational 
maps we were able to create a configuration to study the lattice and perhaps a similar 
technique could be employed to study other crystals. 
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4 Conclusions 



We have introduced an ansatz for a Skyrme lattice with hexagonal symmetry and used 
it as an initial condition in a numerical relaxation of the full nonlinear equations of the 
Skyrme model. The result is further evidence to support the analogy between Skyrmions 
and fullerene shells in carbon chemistry. 

Further work is required to address the issue of the structure of high charge Skyrmions. 
One approach is to collide Skyrmions, but it may also prove useful to use the Skyrme 
lattice, or at least a part of it. The interaction of two such lattices may shed some light 
on the subject of shell formation and it would also be interesting to study the scattering 
of Skyrmions off the lattice. These issues are currently under investigation. 

Given the connection between Skyrmions, monopoles and rational maps it seems 
likely that a similar lattice of monopoles will exist. In this letter we chose a specific elliptic 
function in our ansatz in order to obtain a hexagonal lattice. Other Skyrme fields can be 
obtained which correspond to different lattices by choosing other elliptic functions, though 
we expect them to have higher energy per baryon than the one considered here [|. However, 
in the monopole context all types of lattices would be on an equal footing since monopoles 
are BPS solitons. It would therefore be interesting to see if the Donaldson correspondence 
could be extended to infinite charge monopole lattices, with each characterized by an 
elliptic function. 
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Figure Captions 

Fig. 1. Baryon density isosurface for the Skyrme lattice. The section displayed has baryon 
number four and contains effectively eight hexagons 

Fig. 2. Baryon density isosurface for the Skyrme crystal. Each corner contains a half- 
skyrmion and the total baryon number shown is four. If the threshold for the isosurface 
were to be increased the corners would be connected by links of lower baryon density into 
a crystalline structure. 
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Fig.2. 
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